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Abstract. We study the spectral stability properties of periodic traveling 
waves in the sine-Gordon equation, including viraves of both subluminal and 
supcrluminal propagation velocities as well as waves of both librational and 
rotational types. We prove that only subluminal rotational waves are spectrally 
stable and establish exponential instability in the other three cases. Our proof 
corrects a frequently cited one given by Scott |Sco69b| . 



1. Introduction 

Consider the sine-Gordon equation |Sco69a| in laboratory coordinates, 

utt ~ u^x + sinu ^ 0, (1.1) 

where u is a scalar and {x, t) e M x [0, -l-oo). Although it first appeared in the study 
of the geometry of surfaces with negative Gaussian curvature (see [ETsOOj ) . the sine- 
Gordon equation describes a great variety of physical phenomena as well, such as 
the propagation of magnetic flux on a Josephson line |Sco69aj , elementary particles 
[PS62], modeling of fermions in the Thirring model ColT^, the propagation of 
crystal dislocations |FK39j , and the oscillations of an array of rigid pendula rotating 
under gravity about a common axis with nearest-neighbor torque coupling [Dra83J , 
among others. A comprehensive account of these and other physical applications 
can be found in the review article by Barone et al. |BEMS71| . 

Our paper is concerned with the linearized (spectral) stability of the family of 



periodic traveling wave solutions to the sine-Gordon equation (1.1). To describe 
general traveling wave solutions, one first goes into a frame of reference moving with 
constant velocity c, which in turn amounts to rewriting ( |1.1[ ) in the new independent 
variables z = x — ct and t = t. Thus with u(a;, t) = v{z, t), 

(c^ — l)vzz — '^cVzT + + sin(u) — 0. (1-2) 
In what follows, we will always assume that c 7^ ±1. A traveling wave solution of the 



sine-Gordon equation (1.1 1 is by definition a stationary (r-independent) solution 



of (1.2). Making the ansatz v{z,t) = f{z) implies that / solves the nonlinear 



pendulum equation, that is, the ordinary differential equation 

(c2-l)/"(z)+sin(/(z))=0, (1.3) 
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The pendulum equation ( |1.3[ ) can be integrated once to obtain: 

^(c^-l)/'(z)2 + l-cos(/(z)) = i? (1.4) 

where i? is a constant of integration (the total energy). Clearly, these equations 
are invariant under the shift f{z) i— > f{z) + 27r. 

Definition 1.1 (Wave speed dichotomy of traveling waves). Traveling wave so- 
lutions f{z) with wave speeds satisfying < 1 (respectively > 1) are called 
suhluminal (respectively superluminal) . 

Representative phase portraits corresponding to subluminal and superluminal 
traveling waves are illustrated in Figure [l] It is easy to confirm that except for 



df/dz df/dz 




(a) Subluminal waves: < 1 (b) Superluminal waves: > 1 

Figure 1. Phase portraits of the pendulum equation showing 
both librational waves (closed orbits inside the separatrix) and 
rotational waves (orbits outside the separatrix). The separatrix is 
depicted by a dotted curve. 

solutions corresponding to the separatrices, all traveling wave solutions f{z) are 
periodic modulo 27r. The second dichotomy for such traveling waves is the following. 

Definition 1.2 (Energy dichotomy of periodic traveling waves). Solutions f{z) 



to the pendulum equation (1.3) whose orbits in the phase plane lie outside the 
separatrix are called rotational waves. Solutions whose orbits in the phase plane 
are within the separatrix are called librational waves. 

It is easy to see that librational waves correspond to energies in the range < 
E < 2 and rotational waves correspond to energies with either E < Q (in the 
subluminal case) or E > 2 (in the superluminal case). 

Remark 1.3. In classical mechanics (cf. Goldstein |Gol80| : see also |BM12) ) the 
term libration is borrowed from the astronomical literature describing periodic mo- 
tions in which both the position and the momentum are periodic functions with 
same frequency. The term rotation (sometimes also called circulation or revolu- 
tion) is used to characterize the kind of periodic motion in which the momentum 
is periodic but the position is no longer bounded. Increments by a period in the 
position, however, produce no essential change in the state of the system. In the 
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sequel we will refer to both librational and rotational waves as "periodic traveling 
waves," regardless of whether or not / returns to the same value. 

To study the stability of the periodic traveling wave /, we make the substitution 



w(z,t) :— f{z) + 'w{z,t) in (1.2) and consider the implied behavior of the pertur- 
bation w{z, t). In the simplest approximation as the perturbation is initially small, 
we linearize about ui = to obtain the following: 

(c^ — l)Wzz — 2cWzr + Wtt + COs{f{z))w = 0. (1.5) 

We seek separated solutions with exponential growth rate A G C of the form 

w;(z,r) =p(z)e^^, (1.6) 



which reduces (1.5) to a linear ordinary differential equation for p. 

dz 



p"~2cjXp' + j[X^+cosif{z))]p^0, '-.^^ (P) 



where 7 is defined as: 

7:=^. (1.7) 

— 1 

For all periodic traveling waves f{z), equation ([P]) has periodic coefficients. Let 
T, the fundamental period of /, denote the smallest positive number for which 
/(z + T) =/(z) (mod 2^). 

Definition 1.4. We say that A € C is a temporal eigenvalue if there exists a solution 
to (|P]) which is bounded for all z € M. The set of all temporal eigenvalues is called 
the spectrum of equation ( [P] ), and is denoted fT([P|). 

Definition 1.5. Let f(z) be a periodic traveling wave solution of the sine-Gordon 
equation. We say that f{z) is (temporally spectrally) stable provided that there are 
no temporal eigenvalues in the right half plane. If there is a temporal eigenvalue in 
the right half plane, we say that /(z) is (temporally spectrally) unstable. 

The spectrum it([P]) has the following four- fold symmetry. 

Proposition 1.6. ct([P| = <t^* = -cr^ = -a^* . 

Proof. Suppose A is such that there exists a bounded solution p{z) to equation ( [P| ). 
By taking the complex conjugate of equation (|P]) we see that p(z)* solves equation 
( [P| ) when A = A*. Moreover p{z)* is bounded if p{z) is for all z S M. Thus A* is in 
the spectrum. Likewise for —A, we can consider the function p(— z). Substituting 
these into equation ([P]) gives that —A is in the spectrum as well. Combining these 
last two results completes the proof. □ 



Remark 1.7. The significance of Proposition 1.6 is that the spectral problem (jPj) 
has so-called full Hamiltonian symmetry. Of course this means that spectral insta- 
bility corresponds to the existence of a temporal eigenvalue with nonzero (positive 
or negative) real part. 

Remark 1.8. Under the substitution X = i(^, equation (|P]) can be written in terms 
of a selfadjoint quadratic operator pencil >C(C) given by C{C) ■= C^C2 + C^i + 'Co, 
where 

d d^ 

£2 := -7, £1 := -2c7 • i — , Cq := — + j cos{f{z)). (1.8) 

dz dz'^ 
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Indeed, it is easy to see that (|P]) can be written simply as C{C,)p — 0. In general, a 
spectral problem for a polynomial operator pencil can be reformulated as a genuine 
eigenvalue problem for an operator acting on an appropriate Cartesian product of 
the base space. Upon setting po = p and pi = Xp, we can easily rewrite ([P]) in the 
form 

l.M^aM, Lr.^( V (1.9) 

\PiJ \PiJ V^(c^ - 1)9^^ - cos(/(z)) 2cdzJ 

The matrix operator Lf is frequently called the companion matrix to the pencil 
C. In principle, these observations make available the theory of Krein signatures 
associated to the purely imaginary temporal eigenvalues, which can in turn be 
used to locate non- imaginary points of (t([P|) corresponding to instability. Indeed, 
an alternate proof of Lemma |4.1| below can be given in terms of Krein signature 
theory. See |KM12j and the references therein for further information. 

We are now ready to formulate our main result. 

Theorem 1.9. Subluminal rotational waves are stable. Superluminal rotational 
waves are unstable. Both sub- and superluminal librational waves are unstable. 

In Figure [2] we show numerically computed spectra for each of the four types of 
periodic traveling waves. The structure of (t([P]), including details of the behavior 
near the origin in the A-plane will be described (also for a more general class of 
nonlinear Klein-Gordon type equations) in a subsequent work |JMMP12| . 

Remark 1.10. The complete stability analysis of periodic traveling wave solutions 



of the sine-Gordon equation (1.1) requires two additional nontrivial steps. First, 
one requires some analogue of the spectral theorem guaranteeing completeness of 
the set of nonzero solutions p{z) of equation ([P]) corresponding to A G cr([P]). This 
would allow the general solution t) of the Cauchy initial- value problem for the 
linearized sine-Gordon equation ( |1.5[ ) to be expressed uniquely as an integral over 
(7([P|) of factorized solutions w{z, r) = p{z)e^'^ . Second, one must prove that w{z, r) 



is a good approximation to the solution w(z, r) of (1.2 1 with the same initial data, 
i.e. one must prove nonlinear stability. However, the notion of spectral stability is 
the essential starting point for both of these next steps. 

It should also be said that the kind of stability problem we have in mind from the 
start is the analysis of the effect of perturbations to initial data at r = t = that 
are small in an appropriate function space (e.g. L^(K)). A quite different approach 



would be to exploit the invariance of the sine-Gordon equation (1.1) under the 
Poincare group of Lorentz boosts to make the traveling wave stationary. This 
latter approach has the advantage of avoiding altogether the first-order derivatives 



in equation (|P| arising from the use of a Galilean boost (under which (1.1 1 is 
not invariant) to make the wave stationary; however this method only applies for 
subluminal traveling waves, and moreover physically one is then perturbing the 
traveling wave on a space- like hypersurface different from the line t = 0. 

To our knowledge, Scott |Sco69bj was the first to investigate the stability of 
periodic traveling wave solutions to the sine-Gordon equation. Scott made the 
observation that equation ( [P| ) can be converted by a simple exponential substitution 

9(^)=p(z)e-^^^^ (1.10) 
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Subluminal Rotational 
c = 0.5,£;= -0.5 



Subluminal Librational 
c = 0.5, £; = 0.5 




Superluminal Rotational 
c = 2,£^ = 2.5 
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Figure 2. Numerical plots of the spectrum (t([P]) in the complex 
A-plane representative of the four types of periodic traveling waves. 
Top row: subluminal waves. Bottom row: superluminal waves. 
Left column: rotational waves. Right column: librational waves. 
To compute the part of g([P|) on the imaginary axis it is useful to 
use Corollary |3.2| below to reduce the problem to the calculation 
of the trace of the monodromy matrix for Hill's equation ( |Q| ). To 
compute the non- imaginary part of crW\ one can look for the zero 



level curve of Gp(A) to be defined below (see (3.35)) 



into Hill's equation with spectral parameter /Lt: 



q" + -fCOs{f{z))q = fiq, fi:=j^X'^ 



(Q) 
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This method is advantageous in that Hill's equation has been well-studied, and there 
is a large body of related work dating to the late 19th century (see, for example, 
[MW66| and the references therein). Using this connection and information about 



the spectrum of Hill's equation, Scott formulated the statement of Theorem 1.9 



Unfortunately, the argument in |Sco69bj assumes that the transformation (1.10 1 is 
isospectral, i.e. that the exponential factor e~'^'^^^ is harmless. This is not generally 



true unless A is purely imaginary (as we will prove below; see Corollary 3.2 and 



Lemma 3.3). This led to an erroneous calculation of (t(Ip]) in |Sco69b] . The result 



of Theorem 1.9 is in fact rather well-known in the nonlinear waves community (see, 
for example the influential textbook by Whitham .Whi99] ). however we have not 
been able to find a correct proof in the literature. Our goal in this paper is to give 
a completely rigorous proof of Theorem |1.9[ properly accounting for the effect of 



the exponential factor in the transformation (1.10 1. 

The rest of our paper is organized as follows. In S|2]we define the basic notions 
of spectral analysis for the problems ([P]) and ( [Q| ) . Then in fjs] we present some 
fundamental relations between the spectra of problems ([P]) and ( [Q| ) and we also 
recall the most important properties of Floquet theory for Hill's equation. At this 
point we will have all of the tools in place to give the proof of Theorem |1.9[ all 
details of which are presented in ^|4j Finally, in ^ we give some auxiliary results 
describing aspects of the structure of the spectrum of equatio n ([P] ) in the complex 
A-plane. In particular these results show that the spectrum of (|P[) is not necessarily 
confined to the real and imaginary axes as was suggested in the paper (Sco69b] . 

2. The Eigenvalue Problem 

2.1. Floquet theory for equations ([P]) and ([Q|). Consider a general first-order 
2x2 system with periodic coefficients, of the form 

F'(z;A) = A(z;A)F(z;A), A(z -f T; A) = A(z; A). (2.1) 

Here, the matrix A(z; A) is entire in A for each z e M. The Floquet theory of 
this equation is based on the analysis of the fundamental solution matrix V{z\ A) 



satisfying (2.1) and the initial condition F(0;A) := I. The relevant quantities are 



defined as follows. 

Definition 2.1. The matrix M(A) := F(T; A) is called the monodromy matrix for 
equation ( |2.1[ ) . We denote the trace and determinant of the monodromy matrix by 
A(A) := Mii(A)-hM22(A) and D(A) Mii(A)M22(A)-Mi2(A)M2i(A) respectively. 
The discriminant is the quantity A(A)^ — 4£)(A). 

By Abel's identity D{X) can frequently be calculated in closed form: 

D{X)^exp([ tr(A(z; A))dz ] . (2.2) 



Equation ([Pj) can be written in the form of equation (2.1 ): 

(?)'=^^(^'^)(?)' ^''(^'^)-(-,(A^+cos(/(z))) 2c,x)- (2.3) 
Note that since f{z) is periodic modulo 27r, cos(/(2;)) is also periodic. We denote the 



(2.3) using the subscript "p" . 



quantities defined in Definition 2.1 corresponding to equation (|P|) or equivalently 
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Definition 2.2. The Floquet multipliers p ~ p{\) of equation ([P]) are tlie (spatial) 
eigenvalues of the monodromy matrix Mp(A), i.e., they solve the quadratic equation 
p2- Ap(A)p + i?p(A) = 0. 

The next proposition is a standard result from Floquet theory (e.g. see |Chi99| ) : 

Proposition 2.3. The complex number X is a temporal eigenvalue, i.e. belongs 
to the spectrum (t([P]), if and only if at least one of the Floquet multipliers p has 
modulus I: \p\ — 1. 

For this reason we sometimes refer to as the Floquet spectrum (of ([P])). 

The relation between equations ([P]) and ( [Q| ) given by (1.10) allows us to rein- 



terpret the spectrum criPh as originally defined in Definition 1.4 



Proposition 2.4. The Floquet spectrum g([P|) consists of exactly those A G C such 
that there is a solution q{z) to ([Q|) for which q{z)e^^^^' is bounded for all z G M. 



Just as for ([P]) we can write ( [Q| ) as a system with periodic coefficients: 

^(^=^)-(,-,coV)) J)- ^'-'^ 



We denote the quantities defined in Definition 2.1 corresponding to equation ( [Q| ) 



or equivalently (2.4) using the subscript "g". 

Definition 2.5. The Floquet multipliers 77 = ri[\) of equation ( [Q| ) are the (spatial) 
eigenvalues of the monodromy matrix Mg(A), i.e., they solve the quadratic equation 
7^2- A,(A)?7 + i?,(A) = 0. 



The definition parallel to Definition 1.4 but tailored to ( [Q| ) is the following. 



Definition 2.6. The (Floquet) spectrum g([Q|) of equation ( [Q| ) is the set of A G C 

such that there exists a bounded (on M) solution q{z) to equation ( [Q| ). Equivalently, 
A G cr ([Q| if and only if at least one of the Floquet multipliers -q has modulus 1: 

H = 1. 



From equations (2.3) and (2.4) and formula (2.2) we have that 

Dp(X) = e^"''^^ and ^^(A) = 1. (2.5) 



Remark 2.7. We stress here that in genera l a(|P[) is not the same as g(|Q[). It is an 
elementary fact from Hill's equation theory jMW66) that the spectral problem ( [Q| ) 
is self-adjoint, and hence (7([Q|) is confined to the real and imaginary axes, as /i must 
be real in order for there to exist a bounded solution to equation ( [Q| ). However, it 
is evident from Figure [2] that cr([P]) can have points that are neither purely real nor 
purely imaginary. Moreover, A G (7([P|) corresponds to the existence of a solution 
to ( [Q| ) exhibiting a certain growth rate for z G M, and it is important to keep in 
mind that in determining stability of traveling wave solutions to the sine-Gordon 
equation, it is not the Floquet spectrum of ( [Q| ) per se that will explicitly indicate 
stability, as the latter corresponds to solutions of ([Q]) bounded for z G M. 

On the other hand, the Floquet spectra of ([P]) and ([Q]) are related as we will 
show in 93.11 



Remark 2.8. It is not in general true that the fundamental period T of f will 



coincide with the smallest period of equation (2.3 1. For rotational waves they 



do in fact coincide, but for librational waves it is straightforward to see that T 
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is actually twice the minimal period of cos(/(z)). In computing the monodromy 
matrix M(A) for equations ( [P| ) or ( [Q| , we will always use the fundamental period 
T of /, regardless of whether or not it is also the minimal period of (|P]) or ( [Q| ). 



Remark 2.9. Note that the function Ag(A) as given in Definition 2.1 is frequently 



called the Hill discriminant or sometimes just the discriminant |MW66| . We wish 
to emphasize that as such, this is not the discriminant Ag(A)^ — ADg{X) of the 
characteristic polynomial of the monodromy matrix Mq(A). The reason for our 
terminology is that the discriminant of the characteristic polynomial of Mp(A) plays 
an important role in the analysis that follows but unlike that of Mq(A) it is non- 
trivially related to the corresponding trace, and we wish to use parallel language 
for equations ([P]) and ( [Q| ). 

3. Fundamental Properties of the Floquet Spectra 



The purpose of introducing Hill's equation ( |Q[ ) is that much is known about 
it and its Floquet spectrum. An aim of this section is to state and prove those 
results which give us some information about the structure of the Floquet spectra 
of equations ( [P| ) and ( [Q| ) , but which are not dependent on knowing the type of the 
periodic traveling wave /(z) (i.e. sub- or superluminal, rotational or librational) . 



3.1. Relating the spectra of equations (jPj) and ([Qj). We will begin by relating 
the Floquet multipliers of (|P]) and ( [Q| ). 

Lemma 3.1. The values {p+(A), p_(A)} are the Floquet multipliers of ( [P| ) if and 
only if the corresponding values {r7+(A), ?7_(A)} := {e^'^'''^^p_|_(A), e~'^'''^"^p_ (A)} are 
the Floquet multipliers of ( [Q| ). 

Proof. Let p and q be solutions to (|P]) and ([Q]) respectively, related by equa- 



tion ( |1.10[ ): q = pe-^T^^ We have: 

(|,) = e--H(A) (j) where M(A) (j^^ j) . (3.1) 



Recall the fundamental solution matrices ¥p,q{z; A), solving equations (2.3) or (2.4) 
with the initial conditions Fp ,j(0; A) — I. We let Q{z; A) be the matrix of functions 
given by 

Q(z;A) :=e-^''^^H(A)Fp(z;A). (3.2) 
We observe that Q{z; A) is a matrix solution to ( [Q| ) satisfying Q(0; A) = IHI(A). This 
means that we can write 

Q(z;A) =F,(z;A)H(A). (3.3) 
Equating these last two expressions and evaluating at z = T, and using the fact 
that det(H(A)) = 1, we can then write 

M,(A) = e-"T^'^Il(A)Mp(A)H(A)-^ (3.4) 

The eigenvalues {77+(A), 7?_(A)} of the left-hand side (Floquet multipliers of ([Q])) 
must therefore equal the eigenvalues {e~'^'^'^'^ p+{X), e~'^'^^'^ p-(X)} of the right-hand 
side (where {p+(A), /9_(A)} are the Floquet multipliers of (|P])). □ 

Corollary 3.2. The imaginary number A = i/3, /3 S M, is in the Floquet spectrum 
of ([P|) if and only if it is in the Floquet spectrum of 



The Floquet spectra of (|P]) and ([Q| therefore agree on the imaginary axis. The 
next lemma says that it is the only place where this occurs. 
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Lemma 3.3. Suppose that A G (t^P\ D cr (|Q|). Then A is purely imaginary 



Proof. Suppose that A is in the Floquet spectrum of ( [Q[ ) . Then in particular every 
nonzero solution of ( [Q| ) is bounded away from zero infinitely often as z — > cxd, z S 



Equation (1.10) then implies that every non-zero solution to (^P| is exponentially 



unbounded, and hence A ^ cr([Pj), unless A is purely imaginary. □ 

We now illustrate some of the main differences between g( [P|) a nd (7([Q|). It is 
well-known that if A e cr(|Q|), then both Floquet multipliers of (|Q[) lie on the unit 
circle. However for ([P]) we have instead the following. 

Lemma 3.4. Assume that c ^ (in addition to ^ 1 ). Let A G a([P|) and suppose 
that Re{\) ^ 0. Then one and only one of the Floquet multipliers of ( [P| ) will lie 
on the unit circle. 



Proof. Let p±(A) denote the Floquet multipliers of equation (jPj). Equation (2.5) 
implies that \p+{X)\ ■ \p-{\)\ = 6^^''™'='^). If A e a^, then one of the factors on 
the left hand side is equal to 1. But if Re (A) ^ 0, then the other factor cannot also 
be equal to 1. □ 

Lemma 3.5. Assume that c 7^ (in addition to 7^ 1 ). Suppose A e M\ {0} is in 
the Floquet spectrum of ( [P| ) . Then one of the Floquet multipliers of ( [P| ) satisfies 
— 1, i.e., A is either a periodic (p = or antiperiodic (p — —1) temporal 
eigenvalue. 



Note that by contrast if A e CT([Qj) is real and nonzero it is not in general a periodic 
or antiperiodic Floquet eigenvalue of Hill's equation ([Q|). 



Proof. Equation (2.5) implies that the product of the Floquet multipliers of ([Pj) is 
positive, and is explicitly given by e^'^''^'^ ^ 1. Since A G K, the Floquet multipliers 
are either real or form a complex-conjugate pair. The latter case can be ruled out 
because for A e cr([P]) the product of the (necessarily unimodular) multipliers would 
equal 1. So the multipliers are both real and since A G cr([P]) one of them must be 
±1, in which case the other is given by ■ize^'^''^^ 7^ ±1- D 

We take a moment here to highlight some of the consequences for the Floquet 
spectra of equations (|P]) and ( [Q| ) that are due to the past few results. 

Proposition 3.6. Let \ be a point where the discriminant of ( [P| ) vanishes, i.e. 
the Floquet multipliers are equal. Denote the equal values by p. Then 

(1) The Floquet multipliers of ( [Q| ) are also equal. Denoting the equal values 
by rj, we have that rj = ±1 and correspondingly p = zte'^'^'^'^ . 

(2) A is either real or purely imaginary. Lf it is real and non-zero, then A ^ 
cr([P|), while if it is purely imaginary (including zero), then A G (t([P]). 



Proof. The fact that the Floquet multipliers of ([Qj) are equal follows from Lemma 3.1 
Formula (2.5) implies that 7/^ — Dq{\) — 1, hence rj — ±1. Applying Lemma 3.1 



again implies that p = ±e'^^'^^. Since I77I — 1, we have A G cr(lQ| ), and by self- 
adjointness of ( [Q| ), this means that Im (A^) ~ 0. Corollary |3.2| and Lemma 3.3 
together imply the dichotomy in statement (2). □ 
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3.2. The spectrum of equation ([Qj). A general reference for the Hill's equation 
theory in this section is Magnus and Winkler |MW66| . Consider the Hill's operator 

L:-^+7Cos(/(z)), (3.5) 

where f{z) is the T-periodic (mod 27r) traveling wave solution of the sine-Gordon 
equation about which we are linearizing. The spectrum g(|Q|) is characterized by 
nonzero solutions q of 

Lqiz)^^<^'^q{z), q{z + T) ^ e^U{z) , (3.6) 

where 6* is a real phase angle. (Recall the relation /i = 7^A^; in this section we will 
consider rather than A as the key spectral parameter for Hill's equation ([Q]).) 

The numbers correspond to the periodic eigenvalues of ( [Q| , while the num- 
bers /i'^'^' correspond to the antiperiodic eigenvalues. The set of /i'^-' € M for which 



there exists a nontrivial solution of (3.6) will be denoted E6/(Q). Since cos(/(z)) 



is real, it is clear that I]_6i(Q) = S]£/(Q). The Floquet spectrum cr([Q| is related to 
the union of these sets over 9 as follows: 

7^(10])' = I](Q) y Ee(Q). (3.7) 

— 7r<e<7r 



The ^-Floquet spectrum S(Q) of Hill's operator (3.51 is bounded above. It consists 
of the union of closed intervals 



S(Q) = U [A'2°Vi- A^2:'+2] U [A*^:Vi' A^2n^] (3-' 



»i=0 

where the sequences So(Q) '■= {Mj"^} S7r(Q) :— {Mj'^^} decrease to — oo and 
satisfy the inequalities 

. . . < ^ < ,r < A^f < < ^ < < M^*'^ . (3.9) 
The function /'(z) is a nontrivial T-periodic solution [0 — 0) when /i = 0, that 
is, Lf'{z) = 0, and f'{z + T) = fi^)- Hence one of the periodic eigenvalues ^j"-* 
coincides with /x = 0, and the value of j is determined by oscillation theory (see 
}MW661 Thm. 2.14 (Haupt's Theorem) Q or [CLSSj Thm. 8.3.1]). For librational 
/, f'{z) has exactly two zeros per period and hence either /ij^-* = or /ij"'' = 0. 
On the other hand, for rotational /, f'{z) has no zeros at all and hence /iQ^-* = 0. 
Therefore, we have the following dichotomy for the spectrum of Hill's equation: 

• For librational /, the positive part of the ^-Floquet spectrum E(Q) consists 
of the intervals [^J''l\ U (which may merge into a single 
interval if ^J'^'' = /i2^''). It is possible that — 0, but not necessary; 

otherwise /ij^' = and > 0. 

• For rotational /, the /z-Floquet spectrum S(Q) is a subset of the closed 

negative fi half-line, and /Xg*^^ = belongs to the spectrum. 



"'^The statement of Haupt's theorem in IMW66| contains some typographical errors; the second 
sentence of Theorem 2.14 should be corrected to read "If A = ^'2n—\ oi' - 
exactly 2n — 1 zeros in the half-open interval < a; < tt." 
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In both cases the /i-Floquet spectrum has an unbounded negative part and consists 
of bands separated by gaps. The bands correspond to \Aq\ < 2, and the gaps 



correspond to \Aq\ > 2. Whether or not fj,l 
resolved as foUows. 



(0) 



for hbrational waves can be 



Lemma 3.7. Let f{z) be a Hbrational traveling wave of the sine-Gordon equation. 
Then = > 

Proof. We begin by characterizing the dependence of the fundamental period T 
on c and E. Since f{z) is a Hbrational traveling wave, we have < i? < 2. By 
integrating (1.4) along the appropriate orbits inside the separatrix, we obtain: 

/2(c2 - l)PiE), c2 > 1 

T=< (3.10) 



V2(l-c2)P(2-i?), 



=^<1, 



where 



P{E) 



arccos(l — S) 



df 



(3.11) 



l + cos(/) 

(The subluminal period formula follows from the superluminal period formula by 
the substitution f i-^ f + tt.) 

Now, since A = is a periodic eigenvalue of ( [Q| , we have Aq — 2 at fi — ji^ 
and therefore to prove that ^^2 ^ < 0, it sufhces to show that 

d 



(0) 







< 0. 



(3.12) 



M=0 



In order to establish (|3.12|), we will first show that the monodromy matrix of 
equation 



takes the form: 
M,(0) = 



^Vo{E,c)^{c^-l)T^ 
1 



(3.13) 



where vq{E, c) is a real non-zero constant related to the function f{z) that will be 
defined below, and where Te — dT/dE. 

To prove (3.13), fix a wave speed c ±1 and a total energy E G (0,2). Let 
f{z; E, c) be the periodic wave uniquely determined modulo 2tt from equation ( 1.3 ), 
having total energy E, and for which 

f^{0;E,c)>0 and sin(/(0; c)) = 0. (3.14) 

We claim that when /i = the two-dimensional vector space of solutions to ( |2.4[ ) is 
spanned by the following (subscripts denote partial derivatives): 



and 



fE 
fsz 



(3.15) 



To see this, first differentiate equation (1.3) with respect to either z or E to see 
that both vectors solve equation (2.4) for /i = 0. Then to prove independence, 
differentiate equation (1.4) with respect to E to get 

(c'-l)/./.B = l-/£;sin(/). (3.16) 
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Using this and (1.3 1 gives 



fz fs 

fzz /ez 



fzfsz - fsfzZ = 7 



1 



c2 - 1 



^0. 



We may therefore define a fundamental solution matrix at /i = by: 
Q(z,0) 



fAz-E,c) fE{z;E,c) 
fzz{z;E,c) fEz{z;E,c) 



(3.17) 



(3.18) 



Of course the unique fundamental solution matrix Fg(z, 0) satisfying Fq(0, 0) = I 
can then be expressed as Fq(z,0) = Q(z, 0)Q(0, 0)~^. In particular, by setting 
z = T we can thus write the monodromy matrix at /i = as 

M^{Q) = Fg(r, 0) = Q(T, 0)Q(0, 0)-^. (3.19) 

0). Denote the initial values of 



Next we wish to explicitly relate Q(T, 0) with 
/ and Jz as follows: 



uo{E,c) ■.= f{0;E,c)^f{T;E,c), 
Vo{E,c) := fz{Q;E,c)^ fz{T;E,c). 



(3.20) 



Differentiate the initial conditions (3.20) with respect to E to get that 

OeUo = fsiT; E, c) + TsfziT; E, c) 
^ fEiT;E,c) + TEfME,c) 
= fE{T;E,c)+TEVo{E,c) 



and 



fEz{0;E,c) = dEVo{E,c) ^TE.hAT;E,c) + fEz{T;E,c) 



(3.21) 



(3.22) 



Using periodicity of / and sin(/), equation ( |1.3[ ) and the initial conditions (3.14) 
imply that 



/,,(r;£;,c) = -7sin(/(0;£;,c)) = 0. 
Combining these observations, we conclude: 



Q(o,o) 



fziO;E,c) fE{0;E,c) 
fzziO;E,c) fEz{0;E,c) 

Vo{E,c) dEUo{E,c 
dEVoiE,c 



(3.23) 



(3.24) 



and likewise. 



Q(r,o) = 



fz{T;E,c) fE{T;E,c) 
fzAT;E,c) fEz{T;E,c) 

vq{E,c) dEUo{E,c) -Tevo{E,c) 
dEVo{E,c) 

~TeVo{E,c) 




(3.25) 



= Q(o,o) 



Using these in (3.19) along with det(Q(0,0)) = 7 yields (3.131. 
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Now, differentiatiorj^ of the equation F^(z,A) = Aq{z, X)¥q{z, X) (cf. equation 



(2.4)) with respect to ^ — 7 A gives 



F'gjz, A) - Aq{z, A)F,,4z, A) = A,,^(z, A)F,(z, A) = (^J F,(z, A). (3.26) 

The general solution of this equation is obtained by variation of parameters, that 
is, by the substitution Fg^^(z, A) = Fq(z, A)W(z, A) for some new unknown matrix 
W(z, A). Note that the initial condition Fq(0,A) = L VA implies that W(0, A) = 0. 
It follows easily (also using det(Fq(z, A)) = 1) that 



¥g,^{z,X)=¥g{z,X) 



-F,.n(C,A)F,.i2(C,A) 
]F,4i(C,A)2 



-F,,i2(C,A)2 
F,,n(C,A)F,,i2(C,A) 



(3.27) 



Setting z = T and taking the trace gives 



(M,,22(A)-M,,ii(A)) / F,,n(C,A)F,,i2(C,A)dC 



,2i(A) / F,a2(C,A)^dC + M,^i2(A) / F,,n(C,A)^dC- (3.28) 



Upon taking A = (and hence /i = 0) and using (3.13) we see that 



= -Vo{E,cy{c'~l)TE / F,,n(C,0)^dC 



p.=0 



(3.29) 



Because Fg(2,A) is real-valued for all /i G M and Fq^ii(z,0) does not vanish identi- 
cally for < z < T because Fq,ii(0, 0) = 1, we easily conclude that 



sgn I —A, 



= -sgn((c2-l)r£ 



(3.30) 



d ^ 

where we have used the fact that vq{E, c) 7^ 0. 
Finally, we show that the inequality 

(c^ - 1)Te > (3.31) 

holds strictly for all librational waves regardless of speed. The proof is by direct 
computation of Te. RecaUing (3.10), we compute P'(£') as follows: since < E < 2, 
and hence arccos(l — E) < tt, we can make the substitution w = cos(/) implying 
dw = — sin(/) df = — \/l — w"^ df and hence obtain: 

dw 



^{l - w^){E-l + w)' 
Next, making the substitution w — E{u — 1) + 1, so that dw = E du, yields: 

du 

P(E)^2 , . 

Jo ^-Eiu-l)^-2{u-lWu 

It is now easy to take the derivative with respect to E: 

{u-l)'^du 



P'{E) 



{-E{u-lf -2{u-l)Y2^ 



> 0. 



(3.32) 



(3.33) 



(3.34) 



This immediately implies (3.31 ) and hence completes the proof of the lemma. □ 



■^This part of the argument follows closely |MW66I pp. 15-18]. 
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As a consequence of the fact that ( [Q| ) can be reduced to Lame's equatioij^ the 
structure of the Floquet spectrum S(Q) is even simpler as there is exactly one 
(open) gap. 



(tt) _ (tt) 

^1 Y ^2 ^ Rotational 



(tt) _ (tt) 

^1 ~ ^2 Librational 



Figure 3. A qualitative illustration of the values n for which there 
exists a bounded solution to equation ( [Q|. The fact that there is a 
single gap corresponds to the fact that (|Q[) is reducible to Lame's 
equation. The placement of the band gap in the librational wave 
case is due to the fact that in this case dAq/dfi < at /i = 0. 
The Floquet spectrum of ( [Q| ) in the A-plane is obtained from two 
copies of the /x-spectrum by the relation A = fi^^^ . 

Another consequence of the preceding analysis is that it gives us a qualitative 
picture of the trace Ag for real values of fj,. The following result applies to all four 
types of periodic traveling waves of the sine-Gordon equation: 

Corollary 3.8. For each periodic traveling wave j , there exists a point /i* < for 
which \Aq\ > 2 at fj, = fj,^,, and therefore there exist two nonzero imaginary points 
A = ±z/3*, /?* > 0, /i* = {i'yP.jtY , that are not in the spectrum of ( [Q] ). 

For /X inside the first spectral gap we have that |Aq| > 2. Figure [4] verifies that 
/X = is the largest periodic eigenvalue in the rotational case, while it is the second 
largest in the librational case. 

To finish off this section, we introduce another pair of functions which will prove 
useful in later stability calculations. Let p± (respectively rj±) denote the Floquet 
multipliers of ([P| (respectively of ([Q])). Define the maps Gp_g : C — > M by 

Gp(A):=log|p+(A)|log|p_(A)| 

G,(A):=logh+(A)|logh_(A)|. ^ " > 

Lemma 3.9 (Important properties of Gp and Gq). The functions Gp : C — > M and 
Gq : C — >■ M are continuous, and A G cr([P]) (respectively A £ (7([Q|) ) if and only if 



Indeed, writing p{z) = — ^7 cos(/(2)), we have from l |1.4| that 

{p'(z)f = 4{p(z) - ei){p{z) - e2){p{z) - eg), 

where ei := ^7, 62 := —§7, and 63 := ^^{E — 1), that is, p{z) is a Wcierstrass eUiptic function 
lAkh90' and so we have that ( |Q| can be transformed into Lame's equation (cf. IMW66| and 
[DN74 ). 



PERIODIC TRAVELING WAVES IN THE SINE-GORDON EQUATION 



15 




(a) rotational (subluminal) 



(b) librational (supcrluminal) 



Fig ure 4. The graph of Ag for real values of /i. As in Lemma 
and equation (3.12), it can be seen that /xj"-* — and that 



3J 

d^qjdii < at /i = in the librational wave case. Left panel 
c = 1/^3, E = 4.92. Right panel: c = ^3, = 1 



Gp(A) = (respectively Gq{X) = 0). Also, 

Gq{X) < 0, VA e C 

and 

Gp{X) = (Re (c7AT))2 + GgiX), VA G 



(3.36) 
(3.37) 



Proof. As roots of the characteristic polynomial of a 2 x 2 matrix with entries ana- 
lytic in A, the Floquet multipliers ?7±(A) and p±(A) are analytic with the exception 
of square-root type branch points at the isolated zeros of the associated (entire) dis- 
criminant, and a system of branch cuts connecting them. Moreover, the multipliers 
are well-defined precisely at each such branch point. Upon crossing such a branch 
cut, the two multipliers are simply permuted, and hence any continuous symmet- 
ric function of the two multipliers can be continuously extended to the system of 
branch cuts. It only remains to point out that the logarithms are continuous as the 



multipliers are necessarily nonzero due to (2.5). The statement that the functions 
Gp and Gq detect the Floquet spectra of (|P[) and ( [Q| ) respectively follows directly 
from Proposition |2 . 3 1 an d Definition |2.6[ 

The inequality (3.36) follows from (2.5), which implies that the Floquet multi- 
pliers of ( |Q| ) are reciprocals of each other. The identity (3.37) also follows from the 
reciprocal nature of the Floquet multipliers of (Iq|) together with Lemma 3.1 □ 



We conclude this section by deducing the sign of Gp{X) at certain special points 
in the complex plane. 

Lemma 3.10. For each periodic traveling wave there exists /3* > such that 
Gp{±il3,) < 0. 



Proof. Let > be as in the statement of Corollary 3.8 Since A = i/?* is imagi- 
nary, it follows from ( 3.37[ ) that Gp(i/3*) = Gg(i/?*). But according to Corollary 3.8 
^ o"([Q|) and hence from Lemma [3. 9| we have Gq{i/3^,) ^ 0. But from (3.36) it 
follows that in fact Gq{i/3^) < and hence Gp{ip^) — Gq{if3^) < as desired. □ 
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Lemma 3.11. For each lihrational traveling wave there exists a* € g(|Q[) real and 
positive, and if in addition c 7^ 0, then Gp{a^) > for all such a*. 

Proof. The fact that there exist real and strictly positive points a* in the Floquet 
spectrum of ( |Q| ) for librational waves is a consequence of /ig^'' > and the rela- 
tion ^ = (7A) . To see that Gp{a^) > 0, apply (3.37) with the observation that 
Gq{a^) — and using c ^ 0. □ 

Lemma 3.12. Let Re (A) be sufficiently large in magnitude. Then sgn(Gp(A)) = 
sgn(7), i.e. Gp(A) > for superluminal waves and Gp(A) < for subluminal waves. 

Proof. The coefficient cos(/(z)) is evidently negligible compared with A^ in equation 
([P|), and it can be shown that the Floquet multipliers of (|P]) satisfy 

log(p±(A)) =7AT(c±l) + 0(l), A ^00. (3.38) 



Therefore, directly from the definition (3.351 of Gp(A) we see that 



Gp(A) = [7Re (A)T(c + 1) + 0(1)] [7RC (A)T(c - 1) + 0(1)] 

= 7Re(A)2r2 + 0(Re(A)), A ^ 00, ^ ' ' 

where we have used the definition ( |1.7| of 7. The leading term dominates for Re (A) 
sufficiently large, and therefore the proof is complete. □ 

Remark 3.13. The function Gp is useful not only in studying spectral stability 
of periodic traveling waves, but also in numerically finding the curves of spectrum. 
Indeed Figure [2] was made by finding the zero locus of the function Gp(A) with the 
appropriate parameters. In this way, Gp is a non-analytic but continuous analogue 
of an Evans function capturing in one fell swoop the temporal eigenvalues for all 
real values of the phase 9 of the (unimodular) Floquet multiplier p. This should 
be contrasted with the so-called periodic Evans function F{X;6) |Gar93) . For each 
9 G M., F{X; 9) is an analytic function of A whose (isolated) zeros are those temporal 
eigenvalues for which there exists a unimodular Floquet multiplier p with phase 9. 

4. Stability and Instability of Periodic Traveling Waves 
We are now ready to determine the stability and instability of the four types 



of periodic traveling waves. The proof of Theorem 1.9 follows from Lemma 4.1 



Lemma |4.2[ and Lemma |4.3| that we formulate and prove in this section. 
4.1. Stability. 

Lemma 4.1 (Spectral stability of subluminal rotational waves). Let f{z) be a 
subluminal rotational wave. Then cr(|P]) is purely imaginary. 



Proof. Recall the Hill's operator L defined by (3.5) in [3.2 and the corresponding 



Hill's equation eigenvalue problem (3.6) parametrized by G K. The eigenvalue 
problem (3.6) is selfadjoint for all real 9 and since the resolvent is compact (Green's 
function ge{z,^) — gg{S,,z)* is a continuous and hence Hilbert-Schmidt kernel on 
[0, T]^), it follows from the spectral theorem that for each S M the eigenf unctions 
associated with points p G Se(Q) form an orthonormal basis of L^{0, T). Moreover 
the corresponding generalized Fourier expansion of a smooth function w{z) satisfy- 
ing the condition w{z -\- T) = e^^w{z) is uniformly convergent. From these facts it 
follows that 

{w,Lw) < llwlpmaxEg (Q) VweG(2)(M), w{z + T) = e^ w{z), (4.1) 
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where 



{w, Lw) 



w{z)*Lw{z)dz and 



\w{z)\'^dz > 0. 



It follows from ( |4.1[ ) that for rotational /, 

{w,Lw)<0 VweC(2)(M), w{z + T) = e'^'wiz), 



(4.2) 



(4.3) 



because Se(Q) C M_. That is, L is negative semidefinite in this case (it is strictly 
negative definite for 6' 7^ (mod 27r)). On the other hand, for librational /, L is 
indefinite. 

We re-write the equation (|P|) in terms of Hill's operator L as 



(c2 - l)Lp{z) 



2cA^(z) + A2p(z) =0, 



and thus bounded on M solutions will satisfy 

p{z + T) = e'^p(z), for 



(4.4) 



(4.5) 



For a given £ K, we define (Je(P) C C to be the set of complex A for which 
there exists a nontrivial solution to (4.4) satisfying the boundary condition (4.5 1. 
The Floquet spectrum of ([P]) is the union over 9 of these sets: 

a(P)= y aeiP). 



(4.6) 



-TT<0<7T 



Suppose A e o'(P). Then there exists 9 G M., such that in fact A G a0{P). 
Let p G C^^''(M) denote the corresponding eigenfunction satisfying (4.4| and (4.5). 
Then, multiplying the differential equation through by p{z)* and integrating over 
the fundamental period [0, T] gives 



where 



(c^ - l){p,Lp) - 2imX+\\pfX^ = 0, 



pizrf^{z)dz€ 



(4.7) 



(4.8) 



That m is real follows by integration by parts using the ^-periodicity condition 
satisfied by p. The relation (4.7) can be viewed as a quadratic equation for A. 
Expressing A in terms of {p,Lp), to, and we have: 



A = 



1 



.±^~m^-{c^~l)\\p\mp,Lp) 



(4.9) 



and given the reality of m these values are clearly purely imaginary as long as 
< 1 (the subluminal case) and / is a rotational wave (implying the negative 
semidefiniteness condition {p,Lp) < according to (4.3)). □ 

4.2. Instability. Instability of the sub- and superluminal librational waves and of 
the superluminal rotational waves is, in each case, a consequence of the continuity 
of the function Gp defined by (3.35) and the Intermediate Value Theorem. 

Lemma 4.2 (Spectral instability of librational waves). Let f(z) be a librational 
wave. Then there exists a temporal eigenvalue A* G g([P|) with Rc (A*) > 0. 
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Proof. If c = 0, then it is obvious that (1.10) is the identity transformation and 
hence equations ( [P| and ( [Q| ) coincide implying f7(|P]) = (tQ. Lemma [all] then 
implies the existence of a real positive temporal eigenvalue A* = a* > 0. 

If c 7^ (but 7^ 1), let > be as in the statement of Lemma 3.10 and 
let a* > be as in the statement of Lemma |3.11[ Choose a continuous mapping 
(parametrized curve) A : [0, 1] — C for which A(0) = Aq := «/3* and A(l) = Ai := a*, 
and for which Re(A(t)) > for < t < 1. (To be concrete, one might select the 
straight line X{t) = Ao(l — t) + Ait.) Then Gp{X{t)) is a continuous function 
from [0,1] to M, and Gp(A(0)) < from Lemma [3T0] while G'p(A(l)) > from 
Lemma 13.111 It follows from the Intermediate Value Theorem that there exists 

G (0, 1) such that Gp{X{t^,)) = 0, and therefore A* := A(<*) has a positive real 
part and A* G (t^. □ 

The main idea of the proof is illustrated in the left-hand panel of Figure [5] 




Figure 5. Left: the spectrum o-([Q|) for a librational wave shown 
in red, and a curve A = X{t) connecting Aq — at which Gp < 
with Ai = a* > at which Gp > 0, resulting in a root A* of Gp 
(point of (t(|P])) in the right half-plane. Right: the spectrum it([Q|) 
for a rotational wave shown in red, and a curve A = X{t) connecting 
Ao — j/3* at which Gp < with a point Ai with Re (Ai) 3> 1 at 
which Gp > if the wave is superluminal, resulting in a root A* of 
Gp in the right half-plane. 



Lemma 4.3 (Spectral instability of superluminal rotational waves). Let f{z) be a 
rotational wave and suppose that > 1. Then there exists a temporal eigenvalue 
A* e ct([P]) with Re (A*) > 0. 

Proof. Again take /3* > as in the statement of Lemma [3.10| so that with Aq := 
we have G'p(Ao) < 0. Pick Ai with Re(Ai) > so large that according to 
Lemma 3.12 we have Gp{Xi) > (because > 1). Choosing any continuous 
curve A : [0, 1] ^ C with A(0) = Aq, A(l) = Ai, and Re(A(i)) > for < t < 1 
and applying the Intermediate Value Theorem to Gp{X{t)), the rest of the proof is 
exactly as in that of Lemma |4.2| □ 
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The proof is illustrated in the right-hand panel of Figure [5] Combining Lemma [41] 
Lemma |4.2[ and Lemma |4.3| completes the proof of Theorem |1.9[ 



5. Some Results on the Structure of the Floquet Spectrum 

In this final section, we go beyond the basic (in)stability result of Theorem 1.9 to 
deduce some qualitative features of the Floquet spectrum cr(|P]) that can be easily 
obtained by virtually the same methodology already in place. The first result in 
this direction is the following. 

Proposition 5.1. Let f{z) be a suhluminal librational wave. Then there exists a 
positive real temporal eigenvalue A.^ e cr([P|), and one of the corresponding Floquet 
multipliers satisfies p{X*) = 1 (i.e. A* is a periodic Floquet eigenvalue of ([P|) ). 

Proof. If c = 0, then equations ([P]) and ( [Q| ) coincide and hence the positive real 
periodic eigenvalue A* (Mo^Vt^)^^^ > of ([Q| is also a point of cr(|P]). Therefore 
both Floquet mulipliers of the equivalent problems (|P]) and ( [Q| are equal to 1 for 
A = A* > 0. 

If c ^ 0, let Ao a.ifil^^ ^Y^^ > be as in the statement of Lemma 



so that as Aq is the most positive real point in ct(IQ|) we have Gp(Ao) > 0. Let 



3.11 



Ai > be so large that by Lemma 3.12 we have Gp{\i) < (because the wave is 



subluminal and hence 7 < 0). Then applying the Intermediate Value Theorem to 
the continuous function Gp(Ao(l — t) + Ait) we find a root t^, e (0, 1), and since 
:— Ao(l — t*) + Ai^* > 0, it is a strictly positive real temporal eigenvalue of 



(|Pj). It follows from Lemma 3.5 that A* must be either a periodic or antiperiodic 
temporal eigenvalue of (|P]). To see that A* is in fact a periodic eigenvalue, note 
that since Aq is the most positive real value for which the Floquet multipliers r] 



of Hill's equation ( [Q| ) coincide (with value rj — 1), it follows from Lemma 3.1 



that the Floquet multipliers p of equation (jPj) are coincident at Aq with value 
p = e'^^^oT ^ Q gj^^ g^j-g (Jigtinct for all A > Aq, either forming a distinct real pair or 



a complex-conjugate pair. From (3.38) in the proof of Lemma 3.12 it is clear that 



for sufficiently large real A the multipliers are real and distinct, and it follows that 



they are in fact real and distinct for all A > Aq. From (2.5) it is obvious that the 
product of the real multipliers p is positive for all A S K and therefore they have the 
same definite sign (positive, by consideration of the degenerate value for A = Ao) 
for all A > Aq. Since one of the multipliers p satisfies = 1 for A = A* > Aq it 
follows that in fact p ~ \ > Q. □ 

The next result concerns the nature of the spectrum (t([P]) near the origin in the 
complex A-plane; temporal eigenvalues can only accumulate at A = 0. 

Proposition 5.2. Let f he a librational wave of any speed c ^ ±1, and let g([P|)° 

denote the corresponding Floquet spectrum of ( [P| ) with points on the real and imagi- 
nary X-axes omitted. There is an open neighborhood U of the origin such that A = 
is the only accumulation point o/(t([P|)° n U lying on the real and imaginary axes. 



Proof. By Proposition 1.6 it suffices to consider the part of cr(^P| in the first quad- 



rant of the complex A-planc. Let Aq = i/?* with > being as in the statement of 



Lemma 3.10 and let Ai = a* with a* > being as in the statement of Lemma 3.11 
Let U be an open disk centered at the origin with radius ^ min{|Ao|, |Ai|}. 

We first prove that there exist limit points of cr([P])° on the union of the real and 
imaginary A-axes. Consider the one-parameter family of curves connecting Aq with 
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Ai with parameter p > 0: A = Xp{t) Ao(l - t)P + Xit^ for < i < 1. Exactly 
as in the proof of Lemma |4.2[ the Intermediate Value Theorem provides for each 
p > a point \ — \^,.p G cr( |P[) that is an interior point of the curve A = \p{t) (see 
the left-hand panel of Figure |6|). As a bounded subset of C, the sequence {X*,p}'^=i 
has limit points, and as the curve A = Ap(<) with Q < t < 1 approaches the union 
of the real segment [0, a*] and the imaginary segment [0, as p — >■ cx), the limit 
points of the sequence necessarily lie on this union of segments. 

Now let A denote any limit point of cr([P|)° lying on the union of the segments 
[0,0;*] and [0,i/?*]. Because cr([P]) is a closed set, it follows that A e g([P| ) as well. 
If Im(A) > and Re (A) = 0, then as Im (A) < /3* we have A ^ g-(|Q[) which by 
Corollary 3.2 contradicts the fact that A e (tIIp]). On the other hand, if Re (A) > 



and Im (A) = 0, then as Re (A) < a.^, we have A G g([Q|) which by Lemma 3.3 



contradicts again the fact that A G cr(rP|. It therefore follows that A = 0. □ 




Figure 6. Left: a magnification of the neighborhood of the origin 
from the left-hand panel of Figure [s] showing the family A ~ \it) 
of curves from the proof of Proposition 5.2 and the corresponding 
points A*,p G cr^^ converging to the origin. Right: the local struc- 
ture of the spectrum near the origin for librational waves (both sub- 
and superluminal) consists of a union of two crossing analytic arcs. 



It is an easy consequence of the fact that the elements of the monodromy matrix 
Ip(A) are entire analytic functions of A that for each R > 0, <t^\ H {|A| < R} is a 



finite union of analytic arcs {arcs of spectrum). Together with Proposition |5 . 2| this 
implies that for i? > sufficiently small, (T([P])n{|A| < R} is a finite union of arcs of 
spectrum meeting at the origin and having the full Hamiltonian symmetry implied 



by Proposition 1.6 To deduce further details of the local structure of (t(|P| near 
A = (for example the number of arcs of spectrum meeting at the origin and their 
asymptotic form in the limit of small |A|) requires methods beyond those presented 
in this paper. However, in a forthcoming paper |JMMP12] we carry out a complete 
local analysis by computing the Taylor expansion of M[p(A) about A = 0, and this 
method shows that for librational waves, cr([P]) near the origin consists of exactly 
two curves crossing transversely with opposite slopes as illustrated qualitatively 
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in the right-hand panel of Figure [6j The nature of the Floquet spectrum near the 
origin is of crucial importance in relating spectral stability/instability to Whitham's 
so-called modulational stability criterion |Whi99) (see also [JZlOj ). 
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